Flux fluctuations in a multi-random- walker model and surface growth dynamics 
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We study the dynamics of visitation flux in a multi-random-walker model by comparison to 
surface growth dynamics in which one random walker drops a particle to a node at each time 
the walker visits the node. In each independent experiment (trial or day) for the multi-random- 
walker model, the number of walkers are randomly chosen from the uniform distribution [(Nrw) — 
ANrw, {Nrw} + ANrw]- The averaged fluctuation o(Trw) of the visitations over all nodes i and 
independent experiments is shown to satisfy the power-law dependence on the walk step Trw as 
g{Trw) — Trw 13 - Furthermore two distinct values of the exponent (3 are found on a scale-free 
network, a random network and regular lattices. One is /3j, which is equal to the growth exponent 
P for the surface fluctuation W in one-random- walker model, and the other is /3 = 1. j3i is found 
for small ANrw or for the system governed by the internal intrinsic dynamics. In contrast f3 — 1 is 
found for large ANrw or for the system governed by the external flux variations. The implications 
of our results to the recent studies on fluctuation dynamics of the nodes on networks are discussed. 

PACS numbers: 05.40.-a,89.75.-k,89.75.Da 
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Recent advances of network theories pj have elucidated 
that the dynamical structure of many complex systems 
from the internet trough the biosystems to the social sys- 
tems forms nontrivial networks such as scale-free net- 
works. Therefore the dynamical behaviors of the inter- 
acting (linked) units (nodes) on the networks are very 
important to understand the dynamics of such complex 
systems. Recently an important theory 0, 0, S S on 
the dynamical behaviors on such interacting (correlated) 
nodes is suggested and confirmed on the internet, World 
Wide Web(WWW), river networks and etc. 

The main theory 0, 0] for the dynamics of the node 
is that there exists the power-law relation between the 
average activity or flux (/,;} of each node i to the flux 
fluctuation &; as 



(fi) a - 



(1) 



Furthermore the theory suggested that there can be two 
distinct classes of dynamical systems. One class con- 
sists of the systems with a ~ 1/2 and the other consists 
of those with a — 1. The a oi 1/2 systems are those 
in which the main controlling dynamics is the internal 
intrinsic dynamics. Examples are signal activity in mi- 
croprocessors and time-resolved information in internet 
routers. The a = 1 systems are those in which the main 
dynamics is controlled by the external variations. Exam- 
ples are river networks, highway traffic and World Wide 
Webs. Recently another study Q| suggested a dynam- 
ical model which has nonuniversal a values on complex 
networks by introducing an impact variable to each node. 

To understand the origin of two distinct classes, sev- 
eral dynamical models have been suggested 0, Q • One 
of them is based on a multi-random-walker model. The 
details of the multi-random- walker model are as follows. 
In each experiment or trial d, which was called as "day" 



m Refs. Hi, the number of unbiased random walkers 
Nrw are randomly chosen from the uniform distribution 
[(N RW ) ~ AN RW , {N RW ) + ANrw}- The walkers are 
initially randomly distributed on a given support (net- 
work) with the N nodes (sites). Then the distributed 
walkers simultaneously take preassigned fixed Trw steps 
in a given experiment d and count the number of visits 
(or flux) to the node i during Trw steps. By repeat- 
ing these experiments D times, one can get a series of 
data {fid} with d = 1, 2, D and i = 1, 2, N. From 
the data we can calculate the average flux and the flux 
fluctuation of a node i as 



1 D 

(fi(T RW )) = — fid{T R w) 



d=l 
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erf {Tr W ) = — ^2 fid( T Rw) ~ {fi(T RW )y 



(2) 



(3) 
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^From this model on the scale-free networks with the de- 
gree exponent 7 = 3 the power-law relation <Ti(TRw) — 
{fi{TRw)) a or Eq. (1) was numerically shown. Further- 
more the systems for small ANrw were shown to be- 
long to a ~ 1/2 class, while those for large ANrw were 
shown to satisfy a = 1. ^From these results it was ex- 
plained why a = 1 systems are governed by the external 
flux variations and the behavior for a = 1/2 is from the 
internal intrinsic dynamics 

In the multi-random-walk model, the flux configura- 
tion {fid{TRw)} exactly corresponds to the height con- 
figuration {hid(TRw)} in the following surface growth 
model. In the growth model each of Nrw random walk- 
ers drops a particle on the node i whenever it visits the 
node i and thus fid(TRw) = Ji^Trw)- Recently this 
kind of surface growth model in which Nrw is fixed (not- 
varied) in any experiment d is studied to know the effect 
of random-walk-like colored noises on the dynamical 
scaling properties of the surface roughening y|. In the 
study [3 the support is an ordinary one-dimensional lat- 
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tice and the dynamical scaling of the surface width W 
satisfies the conventional scaling as 



W-L f{T RW /L)-l a z 



(4) 



with z = a s /(3. 

In this paper we want to introduce another method 
to discriminate the systems governed by the internal dy- 
namics from those derived by external input flux vari- 
ations by comparison to corresponding surface growth 
dynamics. The main quantity for the comparison is the 
averaged af(T R w) over the node i, which is defined as 



a 2 (T 



f 



RW) 



N 

N ^ 

i=l 



RW) 



(5) 



As we shall see, (j(Trw) for the multi-random- walker 
model shows the power-law dependence on the walk step 
Trw as a (Trw) — Trw Furthermore, if the governing 
dynamics is not anomalous but the internal intrinsic dy- 
namics, then W(Trw) essentially has the same dynamical 
behavior as that of the mean-square surface fluctuation 
(width) W(T R w) in one-random- walker surface growth 
model. Then Trw is (acts as) the growth time and the 
exponent (3 is (or act as) the growth exponent in dynami- 
cal surface scaling |8(. If W(T R w) in one-random- walker 
model dynamically behaves as W(T RW ) ~ T^ w , then 
&(Trw) for the systems governed by the internal intrin- 
sic dynamics satisfies a (Trw) — ^rw w ith P = Pi- If 
P ^ ft, then governing dynamics should come from the 
external flux variations. 

We first want to show the exact mathematical relation 
of (J to the surface fluctuation W. From Eq. JSJ, 

N 



(J 2 (T RW ) = — y^g-j(TRw-) 



TV ^ 
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Since /y = hid in the corresponding surface growth 
model, the mean-square surface width W 2 is 



W 2 (T RW ) 



1 



d=l 
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(7) 



Thus 



a 2 (T w ) - VF 2 (T W ) = ^E(^E 



(8) 



The difference between a 2 and IF 2 for the given step Trw 
comes from the different orders of taking the average over 
nodes (i.e., 1/N'^2 i fid(TRw)) an d that over experiments 
(or trials) (i.e., l/DJ2 d fid(T RW ))- 

In mormal processes where any particular node does 
not receive extraordinarily smaller or larger flux than 
any other regions, then one can expect that the or- 
der of taking averages makes no anomalous effects and 
thus l/N^.f id (T RW ) ^ l/DJ2 d f ld (T RW ). Then the 
difference a (Trw) — W 2 (Trw) becomes negligible or 
& 2 (Trw) — W 2 (T RW )- Therefore the possible origins 
which make the physically important difference comes 
from the breakdown of the spatial or temporal symme- 
try of the given system. Typical examples for such break- 
downs are the local columnar defect [g, E3 an d the tem- 
porally colored noises H, 0] . When the input flux fluc- 
tuation becomes large or AN R w is large in the multi- 
random-walker model, it is also possible that the input 
flux to some localized region of the support can become 
anomalously larger or smaller and then the difference 
s (Trw) — W 2 (T R w) can have some crucial effects and 
remains for a certain length time. 

To see the explained effects numerically, we analyze 
the numerical data from the simulations for the random- 
walker models. To see the universal features on homoge- 
neous networks as well as inhomogeneous networks, we 
use regular lattices of one and two dimensions as well as 
random network(RN) Q and scale-free network (SFN) 
with 7 = 3 as supports in the simulations. Used size 
or the number TV of nodes (sites) of RN, SFN, and an 
one-dimensional lattice are I0 4 (N = 10 4 ). The used size 
of the two-dimensional square lattice is 100 x 100. Each 
experiment (trial) is repeated over D = 100 times. Ev- 
ery walker in each experiment d takes 10 steps. Data 
for the flux {fid(TRw)} are obtained for the walk steps 
Trw = 1) 2, 10 4 . 

The first simulation data to analyze are those for one- 
random-walkcr model in which only one random walker 
is initially dropped on a randomly chosen node (site) of 
the support in a given experiment. In this model the 
flux dynamics should be governed by the internal intrin- 
sic dynamics and the difference a 2 (Trw) — W 2 (T R w) is 
expected to be negligible as explained in the previous 
paragraph. In Fig. 1 we have plotted the surface fluc- 
tuation W(Trw) an d the average flux variation cf(T R w) 
in one-random-walker model. On the four kinds of sup- 
port (RN, SFN, one and two dimensional lattices) the 
dependence of a (Trw) on Trw (triangles in upper plots 
of Figs. l(a)-(d)) is almost the same as that of W(Trw) 
(straight lines in the same corresponding figures.). We 
can't find any significant difference a 3 (Trw) — W 2 (Trw) ■ 
Furthermore the data for one random walker model sat- 
isfies the scaling ansatz a (Trw) = W(Trw) = T^ w very 
well on both inhomogeneous networks (RN and SFN) and 
homogeneous networks (regular lattices). The estimated 
growth exponents ft for one-random-walker model are 
ft = 0.53(2) for SFN, ft = 0.51(1) for RN, ft = 0.763(1) 
for one dimensional lattices and ft = 0.56(1) for two di- 
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(a) Scale-free network (b) Random network 




(c) 1 d lattice 



(d) 2 d lattice 



FIG. 1: Triangles (a) and straight lines in the upper plots 
of figures (a), (b), (c), and (d) represent the data for the 
average flux fluctuation a (Trw) over nodes and the surface 
fluctuation W(Trw) in the one-random-walker model, respec- 
tively. Other symbols in the upper plots represent 7j(Trw) for 
ANrw = 0,20,40,80,100,200,800,1000,4000,10000 from 
bottom to top in the multi-random walker model. The lower 
plots show the values of exponent j3 for various ANrw ob- 
tained by fitting the data in upper plots to the relation 
&{Trw) — T^ w . The dotted lines in the lower plots de- 
note the values of 0i obtained from the relation W(Trw) — 
~&{Trw) ~ T^ w in one-random-walker model, (a) on Scale- 
free network with 7 = 3. (b) on Random network, (c) on one 
dimensional lattice, (d) on two dimensional lattice. 



mensional lattices. On the inhomogeneous networks such 
as SFN and RN, 0i is very close to = 1/2, which is the 
value of of random ballistic depositions on any support 
. In contrast on one and two dimensional lattices 0i is 
quite different from 1/2. 

In the simulations for multi-random-walker mod- 
els, the number Nrw of walkers in an experiment d 
is selected randomly among the uniform distribution 



[(N RW ) ~ AJVhwt, (N rw ) + AN RW ]. Then the N RW 
walkers are initially randomly distributed over the nodes 
(sites) of the support. In the simulations (Nrw) = 10 4 
is always imposed and ANrw is varied as ANrw — 
0, 20, 40, 80, 100, 200, 800, 10 3 , 4 x 10 3 , 10 4 . Other condi- 
tions are exactly the same as those in the simulations 
for one-random-walker model. The results for the depen- 
dence of a on Trw in multi-random- walker model for var- 
ious ANrw are also shown in upper plots of Figs. 1(a)- 
(d). From the data and the relation a ~ T^ w the ob- 
tained values of for ANrw ' s are displayed in the lower 
plots of Figs. l(a)-(d), where the dotted line denotes the 
values of 0i in one-random-walker model on the corre- 
sponding support. As can be seen from the lower plots, 
value for small ANrw or ANrw < 200 is nearly equal 
to (3i (i.e., = pi) on both homogenous and inhomoge- 
neous networks. On the other hand (3 for large ANrw 
or ANrw > 4 x 10 3 approaches to = 1 (or > 0.9). 
[3 for crossover values of ANrw (or ANrw = 800 or 
1000) has crossover values between {3i and 1. The results 
for the multi-random walker models are as follows. For 
small ANrw or m the systems where the internal intrin- 
sic dynamics is dominant, &(Trw) follows the behavior 
of W(Trw) — Trw m one-random-walker model quite 
well. For large ANrw or m the systems where the exter- 
nal flux fluctuation is dominant for dynamical behaviors, 



a (Trw) follows the behavior cr(T t 



rw) =s Tg w with/? = 1. 

In summary, we have shown that a shows the power- 
law dependence on the walk step Trw as a ~ T^ w in 
multi-random- walker model for various ANrw- To dis- 
criminate the systems governed by the internal intrin- 
sic dynamics from those governed by the external flux 
variations, the dynamical behavior of one-random- walker 
model is suggested as the criterion. If the dependence 
of a in multi-random-walker model on Trw follows the 
dynamical behavior of the surface width in one-random- 
walker model as W(Trw) — T^ w or (3 = 0i, then the 
model is governed by the internal dynamics. This behav- 
ior occurs when ANrw is small. When (3 ^ (3i and 0=1, 
then the model is governed by the external flux varia- 
tions. The — 1 behavior occurs when ANrw is large. 
We also show that this criterion holds for the model not 
only on the inhomogeneous networks but also on homo- 
geneous networks or regular lattices. Furthermore 0i for 
one and two dimensional lattices is greater than 1/2. Es- 
pecially in one dimension the random walker has the tem- 
poral correlation due to the reentrant property and the 
visitation probability to a certain site has some power- 
law correlation in time 0, ■ In two-dimension the reen- 
trant property is marginal [jj. This reentrant property 
of random walker explains 0i > 1/2 in both one and two 
dimensional lattices. In contrast for SFN and RN whose 
effective dimensionality is 00 , the random walkers are 
transparent 9| and the one-random-walker model behave 
as random ballistic deposition and satisfies = 1/2 Q . 

Final comments are on the nonuniversal values of expo- 
nent a for Eq. (1) in regular lattices or homogeneous net- 
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FIG. 2: Dependence of the exponent a on ANrw of multi- 
random- walker model in one (Id) and two (2d) dimensional 
lattices. The data is obtained by fitting a(TRw) versus 
(MTrw)} to the relation (1) or <n(T RW ) ^ (fr{T RW )) a . The 
data for ff%(TRw) an d {fi(TR\v)} in one and two dimensions 
are obtained from the same simulation conditions as those in 
Figs. 1(c) and (d) with a fixed value of Trw as Trw = 100. 



on the degree k(i), the study Q showed nonuniversal a 
values (i.e. 1/2 < a < 1) on scale-free networks. Nonuni- 
versal a values are also possible on the homogeneous 
networks or regular lattices without the impact variable 
V(i). By using the relation (JiiTnw) — {fi{TRw)) a with 
the fixed Trw on one and two dimensional lattices, we 
obtain a values for various A Nrw ■ The results are shown 
in Fig. 2. The simulation conditions except Trw for Fig. 
2 are exactly the same as those for Figs. 1(c) and (d). 
Trw for each experiment (day) is fixed as Trw = 100. 
For ANrw < 500 where the internal intrinsic dynamics 
are dominant, a values are nearly constant. The values 
of the exponent a for small ANrw are a — 0.77(1) in 
one dimension and a = 0.60(2) in two dimension, re- 
spectively. The behavior a > 1/2 for internal intrinsic 
dynamics should come from the reentrant properties of 
random walks. As explained before the reentrant prop- 
erty makes the power-law correlation in time for the prob- 
ability to visit a certain site again. This kind of temporal 
correlation effects makes the probability distribution for 
(fi) deviate from the Gaussian and thus gives nonuni- 
versal a values for small ANrw- For large ANrw & 
approaches 1 even for homogeneous networks as can be 
seen from data for AN RW = 4000, 10000 in Fig. 2. 



works. Recently a study [|| suggested a model in which 
nonuniversal a values (i.e., a ^ 1 or 1/2) for the relation 
Ci(lW) ^ (MTrw)) with the fixed T RW - By using 
the multi-random-walker model with the quenched im- 
pact variable V(i) for each node i, which is dependent 
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